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Abstract 

The Bethe-Salpeter (BS) equation for scalar-scalar bound states in scalar 
theories without derivative coupling is formulated and solved in Minkowski 
space. This is achieved using the perturbation theory integral representation 
(PTIR), which allows these amplitudes to be expressed as integrals over weight 
functions and known singularity structures and hence allows us to convert the 
BS equation into an integral equation involving weight functions. We obtain 
numerical solutions using this formalism for a number of scattering kernels 
to illustrate the generality of the approach. It applies even when the naive 
Wick rotation is invalid. As a check we verify, for example, that this method 
applied to the special case of the massive ladder exchange kernel reproduces 
the same results as are obtained by Wick rotation. 
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I. INTRODUCTION 



There has been considerable recent interest in covariant descriptions of bound states, 
for example, in conjunction with model calculations of high-energy processes such as deep 
inelastic scattering. A fully covariant description of composite bound states is essential for 
the understanding of hadronic structure over the full range of available momentum transfer. 
In a relativistic field theory the two-body component of a bound state is described by the 
appropriate proper (i.e., one-particle irreducible) three-point vertex function or, equivalently, 
by the Bethe-Salpeter (BS) amplitude. An extensive review of the BS equation has been 
given by Nakanishi Other reviews of BS equation studies and many further references 
can be be found in Refs. The BS amplitude satisfies an integral equation whose kernel 

has singularities due to the Minkowski metric. The resultant solutions are not functions 
but mathematical distributions. The singularity structure of these distributions makes it 
difficult to handle the BS equation numerically in Minkowski space. 

We illustrate the BS equation for a scalar theory in Fig. 1, where P) is the BS 
amplitude, P = pi + p2 is the total four-momentum of the bound state and p = rj2Pi — ?7iP2 
is the relative four-momentum for the two scalar constituents. We have then M = \/P^ 
for the bound state mass and also rji + ri2 = I, but otherwise the choice of the two positive 
real numbers rji and 772 is arbitrary. Note that pi = rjiP + p and p2 = TI2P — p. In the 
nonrelativistic limit the natural choice is = ""^1,2/(^1 + ^2), e.g., rji = ri2 = 1/2 for the 
equal mass case. This is the choice that we will make here. 

The renormalized constituent scalar propagators are -D(p^ 2) K{p, q; P) is the renor- 
malized scattering kernel. For example, in simple ladder approximation in a (p'^a model we 
would have K{p, q; P) = {ig){iDa{\p — q]'^)){ig), where D^ij)^) = — + ie) and mo- 

is the (T-particle mass. Note that the corresponding proper (i.e, one-particle irreducible) 
vertex for the bound state is related to the BS amplitude by $ = {iD){iV){iD). We follow 
standard conventions in our definitions of quantities, (see Ref. and also, e.g., [§]). Thus 
the BS equation for any scalar theory can be written as 

^r(pi,P2) = / -^,i^D{ql)) (.r(gi,g2)) mql)) K{p,q;P), (1) 

where similarly to pi and p2 we have defined qi = rjiP + q and q2 = r)2P — q. Equivalently 
in terms of the BS amplitude we can write 

D(pI)-' t(p, P) D(pI)-^ = - / (04 -P) -f^lP. 9! P) (2) 

where the kernel function defined by /(p, g; P) = ~iK{p, q; P) is the form typically used by 
Nakanishi [Q. In ladder approximation for a 0^cr model we see for example that I{p, q; P) = 
/{m?^ — p^ — it). For ease of comparison with earlier work we will use the latter form of 
the BS equation here, i.e., Eq. d^). 

One approach to dealing with the difficulties presented by the Minkowski-space metric is 
to perform an analytic continuation in the relative-energy variable which is the so-called 
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"Wick rotation" ||^. This has been widely used as a means of solving model BS equations, 
e.g., see Ref. |^ and references therein. The special case of the ladder BS equation is solved as 
a function of Euclidean relative momentum in the standard treatment. This is possible since 
in the ladder approximation the kernel is independent of the the total four-momentum P^. 
The difficulties associated with this approach in the general case arise from the fact that since 
the total four-momentum must remain timelike then P ■ p becomes complex. In addition, 
when one uses a "dressed" propagator for the constituent particles or more complicated 
kernels in the BS equation, the validity of the Wick rotation becomes highly nontrivial. 
For example, essentially all of the dressed propagators studied previously in Wick-rotated 
Dyson-Schwinger equation approach contain complex "ghost" poles 0. Hence, the naive 
Wick rotation obtained by the simple transcription of Minkowski metric to the Euclidean 
metric and vice versa is not valid in general. So while the Euclidean-based approaches 
certainly play a very important role and can be useful in model calculations, it is preferable 
to formulate and solve the BS equation directly in Minkowski space. Here we present a 
method to solve the BS equation without Wick rotation by making use of the perturbation 
theory integral representation (PTIR) for the BS amplitude . 

The PTIR method was first used in conjunction with a Wick rotation by Wick and 
Cutkosky for a scalar-scalar bound state with a massless scalar exchange in the ladder ap- 



proximation 0,11]. This is now commonly referred to as the Wick-Cutkosky model. They 



solved the BS equation in terms of a single variable integral representation, which is a spe- 
cial case of the PTIR. Wanders ||12| first introduced the more general two-variable integral 
representation to solve the BS equation for a scalar-scalar bound state with a massive scalar 
exchange in the ladder approximation. It was shown that all invariant scalar-scalar BS am- 
plitudes have this integral representation in the ladder approximation |13;14] and the weight 
function for the corresponding BS amplitude was solved formally by means of Fredholm 
theory |1I^JT5[| . This two- variable integral representation was independently proposed for 
the vertex function on the basis of axiomatic field theory [|T^-|TB|. Following this proposal of 
an integral representation for the vertex function, Nakanishi made a detailed and systematic 



study of the PTIR |10 



The PTIR is a natural extension of the spectral representation for a two-point Green's 
function to an n-point function in a relativistic field theory. Since the Feynman paramet- 
ric integral always exists for any Feynman diagram in perturbative calculations, one can 
always define the integral representation such that the number of independent integration 
parameters is equal to that of invariant squares of external momenta. Each Feynman dia- 
gram contributes to the weight distribution of the parametric integral, so that the complete 
weight function for the renormalized n-point function is identical to the infinite sum of Feyn- 
man diagrams for the renormalized Lagrangian of the theory. Hence, we see that the PTIR of 
a particular renormalized ra-point function is an integral representation of the corresponding 
infinite sum of Feynman diagrams for the renormalized theory with n fixed external lines. 

This paper is organized as follows. In Sec. |I|, we introduce the PTIR of the renormalized 
scattering kernel for scalar models without derivative coupling. We also introduce the PTIR 
for the BS amplitude itself. We discuss the structure of the kernel weight function and 
the BS amplitude. In Sec. ^ an integral equation for the BS amplitude weight function is 
derived and expressed in terms of a new type of kernel function, resulting from the kernel 
weight function and the constituent scalar propagators. The structure of this new kernel 
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function are discussed in Sec. Numerical results are presented in Sec. and we present 
our conclusions and directions for future work in Sec. 



II. PTIR FOR SCALAR THEORIES 

In this treatment we will limit ourselves to studies of bound states of scalar particles 
interacting through a scalar kernel without derivative coupling. Let be the field operator 
for the constituent scalar particles having a renormalized mass m. We will define /i such 
that m + /i is the threshold for particle production in the single channel. For example, 
for a 0^0" model where the renormalized a-particle mass satisfies mo- < 2m or mo- > 2m, 
then we would have n = ma or n = 2m respectively. The renormalized propagator for the 
0-particle can be written in the following spectral form; 



fl(,) = _^J^+r ,a^^), (4) 

\m^ — q^ — ie J{m+ii)^ a — — le ) 

where p^{a) is the renormalized spectral function. Note that P(/,(a) > 0, (see, e.g., 0). 

Following the conventions of Ref. [0], (e.g., pp. 481-487), we can also define centre-of- 
momentum and relative coordinates X = r]iXi+r]2X2 and x = X1 — X2 such that xi = X+r]2X, 
X2 = X — rjix, and P ■ X + p ■ x = pi ■ xi + p2 ■ X2- Hence the Bethe-Salpeter amplitude 
$(p, P) for the bound state of two </)-particles having the total momentum P = Pi + P2 and 
relative momentum p = {ri2Pi — ?7iP2) can be defined as 

(O|T0(xi)0(x2)|P) = e-^^-^(O|T0(r/2a;)0(-r/ix)|P) = e"^^'^ / ^e-^P--^p,P) , (5) 



where we have made use of the translational invariance of the BS amplitude. Equivalently, 
we can write 

$(p,P) = e'^-^y rf'xe'^'-"(O|T0(xi)0(x2)|P) = J d^xe'^-"" {Q\T(f){r]2x)(j){-r]ix)\P) . (6) 

Note that the bound states are normalized such that (P|P') = 2u!p{27t)^6^{P' — P), where 
Up = (P^ + M^)^/^ with M the bound state mass. For a positive energy bound state we 
must have P^ = M^, < P^ < (2m)^, and P° > 0. The normalization condition for the BS 
amplitude is given by 

/^/ 0^,Hq,P)^{D-\pl)D-\pl){2nmp-q) + K{p,q;P)}^^^ 

(7) 

where the conjugate BS amplitude $(j9, P) is defined by 

$(p,P) = e-'^-^y d^xe-'P-^{P\T(j)\xi)(l)\x2)\0) = J d^xe-'P-^{P\T(j)\7]2x)(f)\-7]ix)\0) . 

(8) 
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A. PTIR for Scattering Kernel 



The scattering kernel I{p, q; P) = —iK{p, q; P) describes the process (jxp —>■ 4'4', where p 
and q are the initial and final relative momenta respectively. It is given by the infinite series 
of Feynman diagrams which are two-particle irreducible with respect to the initial and final 
pairs of constituent (p particles. For purely scalar theories without derivative coupling we 
have the formal expression for the full renormalized scattering kernel 



I{p,q;P) = Jd^J di 



n 



+ 



7 



P«s(7,0 



7 



(9) 



where qf is the 4-momentum squared carried by (pi and s,t and u are the usual Mandelstam 
variables. This expression has been derived by Nakanishi (PTIR) |T^. Since only six of these 
are independent due to the relation q\ + ql + q^ + q\ = s + 1 + m, this is also the number of 
independent ^-parameters. Hence, we need only introduce the "mass" parameter 7 and six 
dimensionless Feynman parameters with a constraint. The symbol Vt denotes the integral 
region of such that f2 = {^j | < < 1, = 1 = 1; • • • ? 6)}. The "mass" parameter 7 
represents a spectrum of the scattering kernel. The function Pch(7, gives the weight of the 
spectrum arising from three different channels which can be denoted ch = {st}, {tw}, {us}. 
This PTIR expression follows since any perturbative Feynman diagram for the scattering 
kernel can be written in this form, and hence this must also be true of their sum. In a 
perturbative calculation the weight function pch is calculable as a power series of the cou- 
pling constant for a given interaction Lagrangian. For more general theories involving, e.g., 
fermions and/or derivative couplings, the numerator of Eq. (y) will also contain momenta 
in general. This would considerably complicate the procedure to be presented here and the 
problem of how to extend this approach to those cases is as yet unresolved. One can impose 
additional support properties on the kernel weight functions pc\i{liO by analyzing general 
Feynman integrals for the theory of interest |T^. However, we do not impose any further 



conditions at this stage since we wish to incorporate more general cases such as separable 
kernels, which cannot be written as combinations of ordinary Feynman diagrams. 

To illustrate the approach we will discuss two explicit examples here: 
(a) Scalar-scalar ladder model with massive scalar exchange: The simple t-channel one-a- 
exchange kernel is given by 

i{p, q; p) = 2 ■ (10) 

— [p — q) — 1^ 

The BS equation with this kernel together with perturbative propagator is often referred 
to as the "scalar-scalar ladder model" . The Wick rotated BS equation for this kernel has 
been studied numerically |]T^. We use this kernel as a check of our calculations. 



(b) Generalized ladder kernel: A sum of the one-cr-exchange kernel Eq. ([10|) and a generalized 
kernel with fixed kernel parameter sets {7*-*'', e^*^}. After the Wick rotation this kernel 
becomes complex due to the p ■ P and q ■ P terms, so that solving the BS amplitude as a 
function of Euclidean relative momentum would be very difficult in this case. 
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B. PTIR for BS Amplitude 



Let us now consider the PTIR of the BS amplitude itself. Since the BS amplitude in the 
center-of-momentum rest frame [i.e., with = (M, 0)] forms an irreducible representation 
of the 0(3) rotation group, we can label all of the bound states with the usual 3-dimensional 
angular momentum quantum numbers i and iz- We can thus construct the integral repre- 
sentation of the partial wave BS amplitude in the rest frame and the PTIR automatically 
allows us to then boost to an arbitrary frame. 

For simplicity, let us begin by considering the s-wave amplitude. In order to apply the 
PTIR to the BS amplitude, it is convenient to first consider the equivalent vertex function 
T{p, P) defined such that i^{p, P) = D{{P/2) +p)T{p, P)D{{P/2) -p). The vertex function 
r(j9, P) is a one-particle- irreducible 3-point Green's function, which can be expressed as an 
infinite sum of Feynman diagrams and hence in terms of PTIR has the form derived by 
Nakanishi 
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i i=ii i=i [zi[^ +py + Z2{^ -pY + z^P^\ -le 



In contrast to the case for the scattering kernel, all invariant squares of momenta are in- 
dependent and so a single weight function p^{n,z) is sufficient to describe the sum of all 
allowed FejTiman diagrams. From the perturbative analysis of Feynman graphs for the 
vertex function one obtains the following support property for the vertex weight function 



\2 



P3(K,i') = 0, unless k > max.[{zi + Z2){rn + pY , 

Zi{m + pY + Z2{m - pY + Z3{2mY, 

zi{m - pY + Z2{m + pY + zsi^mY] , (12) 

where "max" is taken separately for each fixed set of z. For the vertex function with a 
timelike total momentum satisfying < P^ < 4nY, it follows that Eq. (|1^) reduces to the 
following two- variable representation 

Tip, P)= [dp [ dC . ^ (13) 



1 ^2 



with the support property 

P3(/3,^) = 0, unless P > (m + p - ^{1 - \C\)VP^ j . (14) 
This two-variable representation was first proposed in axiomatic field theory |17| , [T8 



For simplicity, we will now assume that the constituent 0-propagators have their free 
form, -D(p^) D^ij?) = l/{p^ — m? + ie), or equivalently P(f){a) = in Eq. (^. The 
arguments can be extended in a straightforward way to the case P(f,{a) ^ 0, but we will not 
pursue this possibility here. Then attaching two 0-propagators as external lines yields the 
PTIR for the BS amplitude 
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oo 1 

-i J da J dz 

— oo —1 



V5n(a, z) 



m' 



+ a- (p"^ + zp-P+^) - 



n+2 ' 



(15) 



where we have introduced the non-negative integer dummy parameter n. In order for the 
expression to be meaningful, we see that we must have the boundary conditions 



hm 



i,ra+l 



0, 



. 



(16) 



We see that a partial integration of Eq. ([TsD with respect to a connects weight functions 



with different n, i.e., 



(Pn+i{a,z) = {n + 2) / da'ipn{a\z) 



(17) 



and hence, as stated above, the integer n is an artificial or dummy parameter. While the 
non-negative integer n is arbitrary, a judicious choice can be advantageous in numerical 
calculations, since the larger we take n, the smoother the weight function becomes for a 
given $(p, P). From the support property of p3(/3, z) in Eq. (p!^ , it follows that the weight 
function (^„(a, 2) = when 



a < mm 



0, m + yU 



m + — 
4 



(18) 



This support property can be understood from the relation between p3{P,z) and </)„(«, z). 



I.e., 



00 1 



iPnicy^z) = J d/3 J dC 



PsiPX) 




a ^ 



e R{zX)- 



a 



(19) 



where R{z,z) = [{l — z)/{l—z)]9{z — z) + [{l + z)/{l + z)]9{z — z). Thus the region of support 
for (pn{c(, z), which is specified by the step-functions in Eq. (|19|), is due to the singularities 
of the 0-propagators attached as external lines to the proper bound-state vertex T. Note 
that the integral representation in Eq. (^) and the support property Eq. ([T8|) are valid even 
if we keep the full nonperturbative 0-propagator D{p^) by including the spectral function 
p^{a). In this case Eq. (|T^) for the BS amplitude weight function ipn{o:, z) must be suitably 
generalized and will include an integration over each of the spectral functions from the two 
</)-propagators. 

Now let us consider the BS amplitudes for higher partial waves, i.e, those for bound 
states with non-zero angular momentum (£ > 0). To define these, we consider the rest 
frame of the bound state, i.e., = (M, 0). The momentum-dependent structures, i.e., the 
denominators, of the PTIR's for the bound state proper vertex F and the BS amplitude $ 
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are independent of transformations under the little group belonging to the 3-momentum p 
in this frame. Hence, for higher partial waves with angular momentum quantum number 
£ and third component F and $ can be written as the product of the £-th order solid 
harmonic y^^ip) = IpI^Y^"" (p) and the corresponding PTIR for the scalar (i.e., s-wave) bound 
state in this frame A simple way to understand this is that since the system consists 
only of scalar particles and since in this frame there is only one available 3- vector (i.e., p), 
then higher partial wave amplitudes must be proportional to Yi^(p). The solid harmonics 
are polynomials of their three arguments and since they have the self- reproducing property 
Eq. (B3) we know that in fact the amplitudes must be proportional to y^''{p). Thus 



we have in the rest frame for a bound state with angular momentum quantum numbers 
and 



$t^'^^l(p, P) = -^yHp) [da[ dz- , ^"^"'"^ (20) 

Qth -1 ["^^ + a - {p"^ + zp ■ P + - iej 

where the total momentum P is P = (a/P^, 0) and where we have introduced the lower 
bound ath for the integral range of a, corresponding to the minimal value of the right hand 
side of Eq. (p^. The boundary conditions for (/^[^'(a, z) are identical to those for the s-wave 
Eq. ([161) . The dummy parameter n can always be taken such that the loop-momentum 
integral of the BS equation Eq. (|^) converges. 

By performing a straightforward Lorentz boost we have the integral representation of 
the partial wave BS amplitude in an arbitrary frame 

^ ^ _i [m^ + a~{p^ + zp-P + P2 /4) - ze]"+' ' ^ ^ 

where P is an arbitrary timelike 4- vector with P^ = and p' = A^^{P)p. The Lorentz 
transformation A(P) connects P and the bound-state rest frame 4- vector P' = (M, 0), i.e, 
P = A(P)P'. In the following sections we will study the BS equation Eq. (^ in an arbitrary 
frame in terms of this integral representation. 



III. BS EQUATION FOR THE WEIGHT FUNCTION 

In this section we will reformulate the BS equation Eq. (^ as an integral equation in 
terms of the weight functions. This is the central result of this paper. Substituting the 
integral representation of the partial wave BS amplitude Eq. (pO|) together with the PTIR 
for the scattering kernel Eq. (|^) into the BS equation Eq. (^), the right hand side of the BS 
equation becomes 

oo oo 1 

= -i J2j d-f J J da J d2;pch(7,0v'n 

n -oo -1 
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A 1 

2-k)H 7 - (flch + hchP- q + CchP^ + dch + Cch g ■ -P + /chP ■ -P) - ie 

yJii^ ^ , 



[m"^ + a- {q^ + zq- P + P^/A) - it] 

where {och, fech, Cch, • • • , /ch} are different linear combinations of in each of the three chan- 
nels, ch= {st}, {tu}, {us}. They are listed in Appendix |^. Applying the Feynman param- 
eterization one can perform the g-integral by setting the dummy parameter n such that 
n + 1 > //2. This condition ensures that the g-integral is finite. Due to the self-reproducing 
property of the solid harmonics this integral is again proportional to the solid harmonics 
yl^{Ar^{P)p) (see Appendix 0). Now we multiply the propagators D{P/2 +p)D{P/2 — p) 
into both sides and absorb them into the RHS expression using Feynman parameterization. 
The BS equation then becomes 

oo 1 



$[^A](p^p) = -iy^^''{A~\P)p) J da J dz, 



\ ' [m^ + a-ip^ + zp-P + PV4) - ie]"+^ 

oo 1 

Pch(7,0 J da J dzK^\cx,z-a,z) ip^^{a,z) . (23) 



00 00 1 



n Qth -1 



The kernel function K]^^ {a , z; a , z) is defined with a fixed parameter set of the scattering 
kernel {ach, 6ch, Cch, • • • , /ch} and its definition is given in Appendix y. Comparing Eq. ( |23D 
with Eq. (21), and using the uniqueness theorem of PTIR [|l^, we obtain the following 
integral equation for V5.[f^(Q;, z): 



00 1 



ipl^\a,z) = J da J dz IC^^\a, z;a, z)ip\^\a, z). (24) 



"th -1 



Here we have introduced the total kernel function /C|fl(a, z; a, z), which is the superposition 
of Kl^^{a, z; a, z) with the kernel weight functions Pch(75^) such that 

00 

/CM(a, z;a,z) = Y.JdjJ dip,^{-f, K^^a, z; a, z). (25) 



Q. 

Since the weight functions Pch(7;0 fo^' the scattering kernel are real functions by their 
construction, the total kernel function K!§{a,z]a,z) is real, so that the Eq. ( plj) is a real 
integral equation with two variables a and z. Thus we have transformed the BS equation, 
which is a singular integral equation of complex distributions, into a real integral equation 
which is frame-independent. Once one solves for the BS amplitude weight function, the 
BS amplitude can be written down in an arbitrary frame. This is clearly advantageous for 
applications of the BS amplitude to relativistic problems. As we will see the real weight 
functions can be in fact be real distributions. However, as is evident from Eq. (|T^) if n is 
chosen sufficiently large these can be always transformed into arbitrarily smooth functions 
suitable for numerical solution. 
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The kernel function K^\a, z\ a, z) for a fixed parameter set has the following structure; 

K^n{^^ ^) = ^no5(«)/io'('^5 A + ^a"""^6'(a)/;,|^l(Q;, z) 

—noi^~^\^^{a^ z\ a, z) — gf^i^^i ^! (26) 

The functions z), fcjf' (a, ^; a, z) and g)n^{a, z\ a, z) are defined as Feynman parameter 

integrals in Appendix |C[ The terms containing the function h}^{a^z) are independent of 
z. Since this feature is independent of the kernel parameters, the total kernel function 
/CM(q;, z; a, z) also has this structure for any input scattering kernel. For example, the total 
kernel function with n = can be written as 

/Co'(cS, z\ a, z) = (5(a)7^o'(a, z) — Qo^a, z; a, z), (27) 

where Tig '(a, 2;) and Qq^o, z; a, z) are the kernel-parameter integrals of the functions 
/t|fl(a,z) and gjl'^{a,z;a,z) with the kernel weight function in Eq. (pSj), respectively. This 
structure suggests that the weight function (y9o'(a,z) contains a (5-function 

^^\a,z)= c'5{a)-^^\a,z), (28) 

where c' is a constant and <^|f'(a, 2;) is a function determined by the kernel. If one substitutes 
this form into the PTIR for the BS amplitude Eq. ([T5|), the (5-function term gives the 
product of two propagators with the weight c'/2. As discussed in Section [II B| , this is just a 



consequence of the fact that the BS amplitude contains the kinematical singularity due to 
the two constituent propagators. Thus this part of the structure is totally independent of the 
details of the scattering kernel, (i.e, independent of the scattering kernel weight functions). 

Now, if the function Q^^a, z; a, z) in Eq. (|27| ) vanishes in some region of a around the 
point a = 0, i.e., if the 5-function singularity corresponds to an isolated "pole" contribution 
to the spectrum of the BS amplitude, one can write the homogeneous integral equation 
Eq. (^4D as the following coupled inhomogeneous equations; 

1 00 1 

c J dzn^\a = 0,z) - J da J dzn^\a, z)(f^o\a, z), (29) 



c' 



A 

-1 «th 



1 00 1 



—(fQ\a,z)=c'J dzQQ\a, z; a = 0, z) + J da J dzQQ\a, z;a, z)ipQ\a, z), (30) 
-1 "th -1 

where we have introduced an "eigenvalue" parameter A. Instead of solving for the mass of 
the bound state for a given scattering kernel, it is more convenient to solve the BS equation 
as an "eigenvalue" problem for a fixed bound-state mass parameter P^. We thus solve for 
the eigenvalue A as a function of and the points P^ which give A(P^) = 1 correspond to 
the masses of the bound state. For the weight function (pQ\a,z) containing the 5-function, 
(i.e., the weight function Eq. (^) with non-zero c'), one can rescale the weight function such 
that the strength of the ^-function term is unity. Hence, one can study the bound state 
problem numerically by iteration even for the n = case. Consider the following formal 
expression obtained by iterating the kernel 
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ip§ (a, z) = 5{a) — J dz Qq^ (a, z; a = 0, z) + ■ ■ ■ . (31) 
-1 

There is similarly an expansion for the eigenvalue A 

^1 oo 1 1 

- = J dzn^\a = 0,z) + J da Jdzn^o\c(,z) J dz' g^Q\a, z; a' = 0, z') + ■ ■ ■ , (32) 



-1 



which is the Fredholm series. The scalar-scalar-ladder model has been formally solved by 
Sato by means of the Fredholm solution. He showed that the second iterated kernel is 
bounded and that the Fredholm theorem is applicable for the resulting regular equation 
| T5| . Nakanishi extended his results to arbitrary partial wave solutions [0 . 



To apply rigorous mathematical theorems, such as Fredholm theory, it is necessary to 
know in detail the singularity structure of the kernel function /C|fl(a, z; a, z) for given weight 
functions pch(7;0- Hence, it is very difficult to discuss their applicability for the most 
general form of the scattering kernel. Instead of considering such theorems, we shall take 
the more pragmatic path of establishing whether or not the integral equation Eq. (|2^) is 
numerically tractable. 



IV. STRUCTURE OF THE KERNEL FUNCTION 

In this section we consider the structure of the kernel function given in Eq. (P^j). We first 
consider the possible singularities of the kernel function K]^'^ ( ) for arbitrary n with 

a fixed kernel parameter set (7,0' i-^-, for constant {7, Och, ^ch, Cch, ■ ■ ■ , fch}- We will in this 
section omit for brevity the subscript ch. Apart from the trivial 5-function for n = 0, which 
becomes a step-function for n = 1, the possible singularities of K]^{a^z]a^z) are those of 
the functions /i[f^(a, ^), fc|f^(a, z; a, 2;), and 5'|f^(a, z] a, z). 

Let us first consider the function (^If' (ci, z] a, z). The integration over the Feynman param- 
eter is easily performed due to the ^-function. Reflecting the integral range of the Feynman 
parameter y G [0, 00), we see that g^^{a^ z\ a, z) has finite support as a function of a, a, z 
and z. The function g^^{6i, z\ a, z) becomes after the Feynman parameter integral 



g^^\a,z;a,z) 




2^"+^ ^ e{B,{a, z; a, zf - 4i(a, z- z)C,{cx, z)) 

K=±i {B^{a, z\ a, z)2 - 4y4(a, z\ z)C^{a, z)j 



cz - f + ^ {yt + o) - ^ \ z^ - {e - az) 




where the symbol Pf- stands for the Hadamard finite part. Any singularities that arise 
from the Feynman parameter integrals should be regularized using the Hadamard finite part 
prescription. This is consistent with the ordinary ie prescription for a calculation of Feynman 
diagrams in momentum space |TO,T^. A discussion of the finite parts of singular integrals 
and the detailed calculation of ^J^^re given in Appendix 0. The variables y^ for i = 1,2 
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are the roots of the equation A{a, z; z)y'^ + B^^a, z\ a, z)y + (^^(a, z) = 0. See Appendix 
P for the definition of the functions A{a,z;z), Bf^{a, z; a, z), Cf,{a,z), (f^{z,z) and rj^i^z). 
Since the factor Ck(^, z)y^ + rj^^iz) in the denominator of Eq. (^) is positive definite for any 
positive y^^ this factor does not cause any singularity. Thus for C'k(«, z) < 0, 5'|f^(Q;, z; a, z) 
is regular everywhere except the pole at a = 0, since A{a^ z\ z) is non-negative for any a, 
z and z. For CJ^a^z) > the square root of the following factor in the denominator can 
vanish 

Bi^{a, z; d, zY — AA{a, z\ z)CJypL, z) 
1 



a 



2 

I fl 

2 



a'A{a,z;z)-[^D^{a,z,z)-^C^{a,z)] (34) 



X f aM(a, z] z) - {\lD^{a,z,z) + \fcja^)^ 



2N 



where D^{a,z^z) is the regular function defined in Appendix 0. This singularity occurs 
at the boundary of the step function. Since a square root singularity is integrable, this 
singularity is numerically tractable for Cn{a,z) > 0. In the general case one can apply 
the Hadamard finite part prescription ( P5|) in Appendix for a regularization. However, 
when two roots coincide, namely at the point {a,z) on which Cf^{a,z) = 0, special care 
is necessary to perform the a-integral. One may handle this singularity by the regularized 
expression Eq. ( |D7] ) with the finite regularization parameter e. Thus ip^^^{a,z) may have 
singularities even for a constant kernel parameter set in the general case. Note that this 
singular structure of S'jf' (a, z; a, z) is independent of n and /. 

Now let us turn to the functions h^^{a^z) and k]^{a^z]a^z). Since the step function 
in k]^{a^z]a^z) restricts the integral range of y for a given parameter set, it is enough to 
consider the integral 

In ll/mm, 2/maxj - / ay ^^^^^ ^^^^ ^ ^^^^ ^ ^^^^^ ^6^) 



mm 



in order to discuss the singularity structure. It is easy to show that Eq. (R^) can be written 



as 




mm) i/maxy 



(36) 



with the function 



-^(Z/miri) 2/max) = / dy ^— — — TT/ \ , ■ (^"^J 

' a + y [A{a, z)y^ + B[a, z)y + C\"-+^ 



mm 



The boundary of the integral |/mm and ?/max should be understood as constant, when taking 
the derivative with respect to a and a, although they are fixed as functions of variables by the 
step function of Eq. ( |C16|) . Since the y-integral should be performed using the Hadamard 
finite part prescription, the above operation is valid even though /|^] {ymm, 1/max) contains 
singularities. The integral range [i/min, 1/max] for kj^'^{a,z; a, z) depends on a, z, a and z and 
for z) is fixed to [0, oo). Since A[a, z) > for any z and a that is consistent with the 
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support property of the weight function Eq. (0), the integrand always vanishes sufficiently 
fast for y ^ oo. Thus a singularity occurs only if the set of variables; a, z, a and z, satisfy 
the condition where the denominator of the integrand Ay"^ + By + C vanishes for some y, 
namely B^ - AAC > [[ 

For B^ — AAC > 0, the integration over y yields 



/(l/min, 2/max) = a^A-L+c ji™ ^ (Y^ + t 



2C-~aB , f{l+Xy + e^ 
+ lim — , = m ' 



.0 - AAC \{l-Xy + e\ 



(38) 



with X = "^fj^By^ ^ = (^) i^y"^ + By + C). As shown in Appendix |g the factor 
a^A — aB + C becomes a positive definite function D{z) depending only on z after setting 
a = a. The limit B^ — AAC ^ is also regular, since X is proportional to the factor 
\JB'^ — AAC . Thus a singularity occurs only when 1 ± X or F vanish at the end point of the 
y-integration. Furthermore, neither the derivative with respect to a nor that with respect to 
a generates any new singular point. They simply change the power of the singular behavior. 
Since X is independent of a and ^ = 2^^y, the derivative with respect to a acting on 

-^(z/min, 2/max) docs uot causc strougcr singularities. On the other hand, the derivative with 
respect to a increases the inverse power of Y and 1 ± X. Thus the possible singularity of 
the integral 4^1(2/mm, 2/max) is of the form 

■ (^y^in + By^^ + Cr ^^'(^yLx + ^max + C)-' ^^^^ 

where m is a non-negative integer which does not exceed n and where we have used the fact 
that 1/(1 ± X) oc 1/(1 - X2) oc \l{Ay^ + By + C). It is not difficult to see that the fixed 
end points {y = and y = oo), which are independent of the choice of parameters, cause 
no further singularities provided B{a, z) or C does not vanish. Thus the function /i|f^(a, z) 
has no singularity, if C 7^ 0. On the other hand, fcjfl(a, a, 2;) may have singularities. A 
possible end point of the y-integral, which might give the singularity, is; the positive root 
y^ , {i = 1,2) oi the equation 



A{a,z\z)y + B^{a,z;a,z)y + Cy,{a,z) = (40) 

h zb/2—e+az 

Since the Eq. (1^) can be written as follows 



and the boundary value of y at which k changes sign: ?/bd = + | ^e/-/+b/T *hen ?/bd > 0. 



A( ^ 2 , \ , n CK{z,z)y + rj^{z) ^ 

A{a,z)y + B{a,z)y + C = a, (41) 

1 — 2^ 

and the factor Ck(^, z)y+rjt^{z) is positive definite, the end point singularity aiy = y^ is given 
by Pf ■ l/a". Thus the possible singuralities of /c|f^(«, z; a, z) are Pf ■ l/a"' and Pf ■ l/{Ay1^-\- 



min — Oi ymax) at o — Can be singular. Such a singular point is independent of the choice 
of the variables. One can always avoid it by choosing the dummy parameter to be some n > I + 1. 
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Byhd + C*)"^ with < m < n, where n,m = stands for logarithm. However, kj^{d, z; a, z) 
appears in K]l\a,z;a,z) together with the additional factor na^~^^ so that k]^{a,z;a, z) 
exists only for n > 1. The logarithmic singularities then disappear in K]^{a^ z\ a, z), because 
the derivative with respect to a cancel the logarithmic term in (PSI). Furthermore the factor 
weakens the singularity at a = 0, so that only the singularity Pf ■ 1/a survives for 
n > 1 and any /. 

To summarize, the possible singularities of the kernel function Kl^^(a^z; a, z) for a fixed 
kernel parameter set {och, ^ch, Cch, ■ ■ ■ , fch} are: a trivial 5-function singularity at a = for 
the choice of dummy parameter n = 0; Pf ■ l/ct for tt, > 1; Pf ■ l/{Ay^^ + -B?/bd + C*)™ with 
1 < m < n for n > 1; the singularity of Eq. (|3^ at the boundary of the support for all n 
and /. 



V. NUMERICAL RESULTS 

In this section we present numerical solutions for the BS amplitude for bound states in 



scalar theories using Eq.(p4D for two simple choices of scattering kernels: (a) pure ladder 
kernel with massive scalar exchange, and (b) a generalized kernel combined with the above 
pure ladder kernel. 

The scattering kernel (a), i.e., the one-u-exchange kernel depending only on t = {p — qY, 
is given by Eq. ([T0|) . This corresponds to choosing for the kernel in Eq. (H) say: ptu = Pus = 
and in the st-channel 7 = ml, and ast = Cst = 1, bst = -2, 4t = est = fst = 0, [c.f., Eq. (||)]. 
This corresponds to choosing pst to be some appropriate product of 5-functions multiplied 
by g"^. In the pure ladder case it is convenient (and traditional) to factorize out the coupling 
constant g'^ and a factor of (47r)^, by defining the "eigenvalue" A = (7^/(47r)^ []T3|,|T^,|Ty[ . 
Thus it is usual to fix the bound state mass and then to solve for the coupling g"^, which 
is what we have done here. In the numerical calculations it is necessary to regularize the 
integrable square root singularity discussed in the previous section (see Eq. (|3^), i.e., we 
kept the regularization parameter e in the Hadamard finite part Pf-l/x" small but finite, 
typically ~ 10~^ or less. Of course, we verified that our solutions were independent of e 
provided it was chosen sufficiently small. We began by discretizing the a and z axes and 
then solved the integral equation by iteration from some initial assumed weight function. It 
was also, of course, confirmed that the solution was robustly independent of the choice of 
starting weight function, the number of grid points, and the maximum grid value of a (when 
the latter two of these were chosen suitably large). We have solved over a range of and 
TTLrj and all solutions reproduced the eigenvalues obtained in Euclidean space after a Wick 
rotation by Linden and Mitter |^ . We found that this simple method produces convergence 
in typically 10 cycles or so and for the relatively modest grid choices used gave results which 
agreed to within approximately 1% of the Wick-rotated ones. In Fig. (|^) we plot the weight 
function 0„=i(q;,2;) for the bound state with mass given by P'^/m? = {M/mY = 3.24, 
where m is the mass of the constituent particles and where the mass of the exchanged 
particle is rrifj/m = 0.5. This solution was generated by solving for n = and then using 
Eq. ( [T7| ) to convert to the n = 1 form of the solution. This has the advantage of filtering 
out numerical noise and the resulting solution is more readily graphically represented. The 
coupling constant for this case was found to be g'^/lAnY = 1.02, c.f., the Linden and Mitter 
result of 1.0349. 
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For a nontrivial example of the scattering kernel, where the naive Wick-rotation is not 
possible, we studied the following simple generalization of the pure ladder case: a sum of the 
one-(T-exchange kernel (a) and the generalized kernel term with the following fixed kernel 
parameter set 

a,t = 0.47261150181, 
6st = -0.29743163287, 
c,t = 0.58277042955, 
4t = 0.28282145969, 
est = -0.23965580016, 
/st = 0.32196629047, 

7/m2 = 2.25, (42) 

weighted with a factor (l/4)[5f^/(47r)^]. In addition, the resulting total kernel was sym- 
metrized in ^ ^ —z so that we obtain appropriate symmetry for s-wave normal solutions. 
We obtained these parameters by generating in the st-channel with a random number 
generator. We then solved for A = (7^/(47r)^ as for the pure ladder case but for this case the 
n = 1 was solved for directly, without first obtaining the n = solution. In Fig. (|^) we plot 
the weight function 0„=i(a;, z) with the mass parameters again chosen as P'^/m'^ = 3.24 and 
m^/m = 0.5. Some additional structure is apparent, as expected, and some residual numer- 
ical noise is also present in this figure. Methods of reducing this numerical noise by finding 
more sophisticated choices of grid and integration method are being investigated ||20|. The 



resulting eigenvalue is A = 0.889, which shows that the addition to the pure ladder kernel 
has enhanced the binding and so is attractive. We found that the convergence and the 
stability of the eigenvalues with varying the number of grid points, e, and a^^x are excellent 
as for the pure ladder case above. 



VI. SUMMARY AND CONCLUSIONS 

We have derived a real integral equation for the weight function of the scalar-scalar Bethe- 
Salpeter (BS) amplitude from the BS equation for scalar theories without derivative coupling. 
This was achieved using the perturbation theory integral representation (PTIR), which is 
an extension of the spectral representation for two-point Green's functions, for both the 
scattering kernel [Eq. (^] and the BS amplitude itself [Eqs. (^,^)]. The uniqueness theorem 
of the PTIR and the appropriate application of Feynman parameterization then led to the 
central result of the paper given in Eq. (P^. We demonstrated that this integral equation is 
numerically tractable for both the pure ladder case and an arbitrary generalization of this. 
We have verified that our numerical solutions agree with those previously obtained from a 
Euclidean treatment of the pure ladder limit. 

This represents a potentially powerful new approach to obtaining solutions of the BS 



equation and additional results and applications are currently being investigated |^. The 
separable kernel case should be studied, since it has exact solutions and so is a further inde- 
pendent test of the approach developed here. While our detailed discussions were limited to 
the equal mass case, it is worthwhile to investigate generalizations including nonperturbative 
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constituent propagators and the heavy-light bound state hmit to see under which conditions 
an approximate Klein-Gordon equation can result. It is also important to find a means to 
generalize this approach to include derivative couplings and fermions. 
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APPENDIX A: PTIR FOR SCATTERING KERNEL 

In this appendix we list the dimensionless coefficients {act, Kh, Cch, • • • , fch} in Eq. (p2D 
for different channels ch= {st}, {tu}, {us} in terms of the Feynman parameters defined in 
Eq. (D. 



st tu us 





6 + 6 + 6 


6 + 6 + 6 + 6 


6 + 6 + 6 


bch 


-26 


2(6-6) 


26 


Cch 


6 + 6 + 6 


6 + 6 + 6 + 6 


6 + 6 + 6 


dch 


1(6 + 6 + 6 + 6) +6 


i(6 + 6 + 6 + 6) 


i(6 + 6 + 6 + 6) + e 


Cch 


6-6 


6-6 


6-6 


fch 


6-6 


6-6 


6-6 



APPENDIX B: SOLID HARMONICS 



Let us consider the integral 



I{n,mJ;p,P) = J 



{2tt)H [7 - (ag2 j^bp - q + cp'^ + dP'^ + eq- P + fp- P)- ie] 
ye'{A-\P)q) 



^[m^ + a-{q'^ + zq-P + P^/A) - ie] 
Applying the Feynman parameterization we have 

/(n, m, £; p, P) 

d'^q T(n + m) } , x'^'Hl - x)""'^ 
' ax- 



(Bl) 



{27T)HT{n)T{m) (x + (1 - x)a)"+^ 



(B2) 



yl^iA'\P)q) 



(2,'+(l~x)a)2-'-' 



H 



+ m^-(p'^ + ^ + Jp-P 



(xz+{l-x)e)P+(l-x)bp 
2{x+{l~x)a) 



n+m 1 



where F, H and J are 
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F = z)-^ + B{a, z) + C^^, 
1 — X X 

X 

J = f-b/2z + ^{af-el), 
X 2 

and A{a,z), B{a,z) and C are defined in Appendix |y. Introducing a new variable 

q' = A-\P) (q + ^"'+^2(x+a-.?a)""^'"' ) and recalling the fact A-\P)P = {Vl^,0), one can 
factorize the solid harmonics outside the integral by making use of the following property of 
the solid harmonics; 

J d'qF{f)yHq + P) = yhp) J d'qFif), (B3) 

where F is a sufficiently rapidly decreasing function which gives a finite integral [|T^]. The 
integration over q then yields 

J(n, m, £; p, P) 



r(ra)r(m) (47r)2 \ 2 / 7 [cy + A]"+™ 

^ ^ A{a,z)y^+B{a,z)y+C , _2 ' in+m-2- (^4) 



le 



where we have introduced a new integration variable y such that y = x/ (1— x). Thus the inte- 
gration over the loop momentum reproduces the same solid harmonic. This self-reproducing 
property of the solid harmonics ensures that the integral representation, Eq. (|2l|), forms an 
irreducible representation of the Poincare group. 



APPENDIX C: KERNEL FUNCTION 



In this Appendix we present the explicit expression of the kernel function for the constant 
kernel parameter set a, 6, c, d, e, /. The kernel function Kn{a, z; a, z) for the weight function 
Lp]^{a^z) is defined as the following Feynman parameter integral 



K^^\a,z]a,z) 



1 1 



(47r)2 2 
oa 



dy 



y 



n+l. 



a + y) 



[A{a,z)y' + B{a,z)y + Cr^' 

A{a, z)y^ + B{a, z)y + C' 



(CI) 



e{a)-e{a- R{z,G{z)) 



cy + A 

where have suppressed the explicit dependence on the kernel parameters in A{a, z),B{a, z) 
and C. These are defined as 



z) 



a + m 



p2 



B{a, z) = aa + 'J + {a 



P2 



m 



{Ad-2ezy- 



(C2) 



C = a-f-A[m' - — ] - {Aad - e^)^ 
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where A = ac — h"^ /A. The functions R{z^ z) and G{z; y) are defined as 



R(z.z) = !^e(z-z) + ^e(z-z), (C3) 

Note that these functions are bounded, such that 

0<R{z,z)<l for |^|,|z|<l, (C5) 
\G{z;y)\<l for |z| < 1 and y > 0, (C6) 

provided \ f ± b/2\ < c and \af — eb/2\ < A. These conditions are automatically satisfied, 
which is readily seen if one rewrites a, 6, . . . , / by ^ G fi, with ^ = {^i\0 < < 1, = 
1 (i = 1, . . . ,6)}. To integrate over y, the denominators of the integrand [i.e., y"-~^''- and 
{A{a,z)y'^ + B{a,z)y + C)~^^~^^^] should be understood as the Hadamard finite part, if 
necessary. Since G{z = 0;y) = 0, the kernel function Kjl^^a, z; a, z) vanish identically at 
z = 0, so that we need only consider the case z from now on. 

Performing the derivative with respect to a, the kernel function Eq. (|Cl| ) becomes 

K^\a,z;a,z) = ^{a-e{a))h^^\a,z) 



na^ ^k]^\a,z;a,z) + g^^\a,z;a,z) (C7) 



with the functions z), A;jfl(a, z; a, z) and z] a, z); 



1 


-1 


(47r; 




1 


-i 


(47r; 


22 


xe 1 


a — 


1 


a" 


(4vr; 


2 2 


x6 I 


a — 



fey y^'+^a + yy 

2) Jo ^[A{a,z)y' + B{a,z)y + Cr 



+ A / 

mi- - ^ X - , feV y-+\a + yr-^-' 

q^Ha, z; a, z)= - — -- — — / ay — rr 

' ' ^ ^ ' 2^ io ^[A{a,z)y'' + B{a,z)y + Cr^^ 

A{a,z)y^ + B{a,z) + C\ 

z, G{z; y)) — . 

cy + A / 



(C9) 



(CIO) 



Now let us consider the ?/-integration. We first rewrite the step function and the delta 
function as functions of y. Substituting ( P3|) and (P^, the argument of the step function 
and the (5-function can be written as 

- A{a,z)y'^ + B{a,z)y + C 
a - R{z,G{z;y)) 



cy + A 

A{a, z; z)y^ + B^{a, z; a, z)y + (^^(a, z) 
CK{z,z)y + ri^{z) 



(Cll) 
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with K = ±1 for 



cz-f + ^z\{y + a)-^ l^z ~ {e - azU J 



(C12) 



and K = when LHS of ( |U12D vanishes. We have introduced the coefficients for the denom- 
inator defined as 

Uz, z)=(^c + K (^^z -/))(! + '^^) + (1 - - /) ^' 

Vn{z) =(^A-K(^af- ^e^ ^ (1 + kz) + (1 - k') (^^e - af^ z, (C13) 

and for the numerator we have 

A{oL, z; ^) = (1 — z^)y4(a, z) , 

z] a, ^) = (1 - z^)B{a, z) - Ck{z, z)a , 

C^{a,z) = {l-z^)C-r,^{-z)a. (CM) 

Note that C,k,{z,z) > and ^ for \z\, \z\ < 1. Thus the denominator is non-negative 

for positive y, so that the sign of the argument ( |C11| ) becomes positive only if the following 
equation has a positive root 

A{a, z; z)y'^ + B^{a, z; a, z)y + C^{a, ^) = , (C15) 

which requires that )2 - 4A{a,z;z)C^{a,z) > 0. Thus the step function in 

) becomes 

A{a,z)y'^ + B{a,z) + C\ 



e {a-R{z,G{z; y))'- 



cy + A 



e{B^{a, z; a, zf - AA{a, z; z)) e{y - y-^)e{y^ - y) 



K=±l 



xe{K [{cz - / + b/2){y + a) - 6/2 {zb/2 - (e - az))]) 
where y^ with i = 1,2 are two roots of the Eq. (|C15|) 



2/2 



B^,{a,z;a,z) T \^ B^{a, z; a, z^ -4:A{a,z; z)C^{a,z) 



(C16) 



(C17) 



2A{a, z; z) 

The ^/-integration for the functions h^^{a,z), k^^{a,z;a,z) and g^^{a,z;a, z) are discussed 
in Appendix 0. 

Finally we list some useful linear combinations of the functions A{a, z) and B{a, z); 



D{z) = a^A{a,z) - aB{a,z) + C 

= T(^^ -— J+(e-a^) — , 

Dfi{a, z, z)= o?A{a, z\ z) — aB^{a, z; a, z) + 6*^(0;, z) 
= (1 - z^)D{z) + (aC(^, z) - ri^{z)) a 



(C18) 



(1 - z'')D{z) + 



-(1 + nz) — (e — az){z + k) 



a 



(C19) 
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APPENDIX D: FINITE PART OF SINGULAR INTEGRALS 



In this Appendix we review the Hadamard finite part prescription and present the deriva- 
tion of Eqs. ( p3D and (pSl). We use the following definition of the Hadamard finite part 
Pf ■ 1/a;" together with the n-th derivative of the 5-function 



lim ■ 



1 



Pf ■ l/x'^TiT^- 



-1 



in— 1 



for n = 1,2, 



e^o(a;±ze)" ' (n-1)! 

The following formulae are useful; 
d 



X) 



, Pf ■ = (-n)Pf ■ 1/x 
ax 

x6'{x) = —S{x), 

x'^pf ■ = Pf ■ 



n+l 



(Dl) 



(D2) 

(D3) 
(D4) 



for a positive integer k. These are a generalization of original definition of the Hadamard 
finite part 



//■OO 
dx Pf ■ x^9{x)f{x) = lim J x^f{x) + ^ 



(D5) 



where ReA + A; + 2>0 and A is not a negative integer. We also apply this prescription for 
the function g]^{a, z; a, z), as necessary. 

Having understood the (5-function as the limit of Eq. (pi|) , we shall now derive the 
Eq. ( p3D from Eq. ( |C10| ). Recalling the Eq. the delta function can be written as 



5 — R{z, G{z; y)) 
C^{z,z)y + r]^{z 



A{a,z)y^ + B{a,z)y + C' 



cy + A 



Im— 



A{a, z; z)y'^ + B^{a, z; a, z)y + C^,{a, z) - ie 



(D6) 



where we have used the fact Ck(-2) z)y + rii^{z) > 0. Thus the function gi^^ia, z; a, z) is written 
as the e — limit 



g'^^\a,z;a,z) 



(47r 



^ pf.i 



a 



:i--z 



2\n+l 



^ ^.^ e{B^{a, z- a, -zf - AA{a, z; z)C^{a, z)) ^ yr^' ^g + y^y-i-i 
K=±i (B^{a, z- a, zf - 4i(a, z; z)[C^[ci, z) - ze)) ^=1,2 Avt + ^k(^))" 



- / + - (yf + a) - - z 



[e — az] 



where y^ are root of the equation; 

A{pL, z] z)y'^ + z] a, z)y + (^^(a, z) — ie = 0. 

It is now easy to rewrite this expression of the form Eq. 



(D7) 



(D8) 
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Now, let us turn to the evaluation of the integral Eq. 



mm 



Changing the integration variable such that 



^=\^] + % + (DIO) 

2C + By ' ^ ' 



the integral for B — AAC > can be written as follows; 

-?^(l/min, 1/max) = a'^A-dB+C | 4 Jy dY Pi ■ — 



^max 



+ '^-'^ r rfXfPf.— Pf.^)L(D12) 

where Y^i^, X^i^i, ■ ■ ■ are the values of Y and X a.t y = ymin, Umax- It is easy to derive the 
Eq. (|38|) by performing the integral over X and Y with the definition of the Hadamard finite 
part ( pri) . 
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FIGURES 



FIG. 1. The Bethe-Salpeter (BS) equation for two scalar constituents in terms of the bound 
state proper vertex (a) and in terms of the BS amplitude (b). 



FIG. 2. The s-wave {i = 0) Bethe-Salpeter (BS) amplitude weight function for n = 1, 
(f)i~^\a, z), for the pure ladder kernel function described in the text. 



FIG. 3. The s-wave {i = 0) Bethe-Salpeter (BS) amplitude weight function for n = 1, 
0f~°'(a, z), for the example general kernel function described in the text. 
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